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We study vison excitations in a quantum dimer model interpolating between the Rokhsar-Kivelson 
models on the square and triangular lattices. In the square-lattice case, the model is known to be 
critical and characterized by U(l) topological quantum numbers. Introducing diagonal dimers brings 
the model to a Z2 resonating-valence-bond phase. We study variationally the emergence of vison 
excitations at low concentration of diagonal dimers, close to the critical point. We find that, in 
this regime, vison excitations are large in size and their structure resembles vortices in type-II 
superconductors. 



I. INTRODUCTION 

The resonating-valence-bond (RVB) state is one of the 
most exciting proposals for strongly correlated phases in 
two dimensions. Of potential relevance to frustrate d sp in 
systems^' 2 -' anc ] high-temperature superconductivity^, it 
has evaded a direct identification in experiments and in 
studies of realistic models. However, a spinless analogue 
of the RVB phase can be systematically studied in quan- 
tum dimer models, where this phase is rigorously con- 
firmed and a cces sible to a variety of analytic and numer- 
ical methods^'. Furthermore, quantum dimer models 
have also motivated the construction of SU(2) invariant 
examples of RVB phases in spinful sytstems-^. 

The main strongly correlated feature of the RVB phase 
is the existence of fractionalized excitations. In a spin 
system, in the conjectured RVB phase, two types of ele- 
mentary excitations must be present: spinons (spin-1/2 
excitations) and visons (Z2 vortices). They have rela- 
tive semionic statistics: a —1 fact or for an elementary 
braiding of a spinon around a visonLLiiLUJ i n dimer mod- 
els, spinon excitations are absent (or effectively pushed 
to infinitely high energy)- 1 , and only vison excitations 
and their combinations appear in the spectrum^. An- 
other consequence of the existence of vison excitations 
is the so called "topological order" : a degeneracy of the 
ground states (and low-lying excitations) for systems on 
multiply connected domain^ 12 * 15 !. 

An interesting question arises regarding the fate of 
these elementary excitations near phase transitions from 
the RVB phase to neighboring non-topological phases. 
Two generic phase transitions of this sort may be envi- 
sioned. Either singlets (or dimers in dimer models) or- 
der into a valence-bond-crystal, or spins order into some 
sort of magnetic order. These two possible paths for the 
RVB state were identified in Ref.Ql)]as two conditions for 
the topological order in Gutzwiller-projected wave func- 
tions. A field-theoretic model of these phase transitions 
has been proposed and analyzed in Ref. [T7J 

In the present paper, we consider another, non-generic 
phase transition into a U(l) critical state. This transi- 



tion is most easily realized by deforming a model on a 
bipartite lattice (with U(l) winding numbers) into one 
on a non-bipartite lattice (with the U(l) symmetry bro- 
ken down to Z2). While this transition represents a non- 
generic, "fine-tuned" situation, it also allows for an an- 
alytic treatment of quantum dynamics near the critical 
point. Specifically, we consider an interpolation between 
the Rokhsar-Kivelson (RK) models on the square and 
triangular lattices in such a way that the exact solvabil- 
ity of the RK ground states is preserved throughout the 
interpolation. 

In such a model, close to the critical point, we can vari- 
ationally construct vison excitations, thus estimating the 
size of the vison gap and identifying the associated length 
scale (the vison size). We find that the vison gap closes 
at the transition point while the vison size diverges. The 
shape of the vison in this limit resembles vortices in type- 
II superconductors, with the vison size corresponding to 
the magnetic penetration length. 

The paper is organized as follows. In Section II, we 
briefly review the properties of the RK dimer models on 
the square and triangular lattices. In Section III, we 
describe the interpolating model which connects those 
lattices. In Section IV, we construct variational vison 
excitations near the critical point and discuss their prop- 
erties. Finally, in Section V we discuss some interesting 
aspects of our results. 

II. ROKHSAR-KIVELSON DIMER MODELS 
ON SQUARE AND TRIANGULAR LATTICES 

The Rokhsar-Kivelson dimer model on any lattice is 
defined in the Hilbert space of all fully pac ked c overings 
of the lattice by dimers by the Hamiltoniar i 12 * 18 ! 

^RK-E(- < l«)(Hl+H"X"l) (1) 

where the sum is taken over all tetragonal plaquettes of 
the lattice. In the present work, we are interested in 
this model on the square and triangular lattices in two 
dimensions, in which cases the sums are understood as 



those over all the plaquettes of the square lattice and over 
all the two-plaquette rhombi of the triangular lattice, re- 
spectively. 

On both lattices, there is a special case of the model 
([I]) with v = t (called the RK point), in which the ground 
state is exactly known and given by all the allowed dimer 
configurations with equal amplitudes^, 

l GS > = TfEl c > ( 2 ) 

c 

(here the sum is taken over all dimer configurations, and 
Z is a normalization factor). However, the properties of 
this state and of the lowest excitations above it are very 
different in the cases of square and triangular lattices, 
due to the fact that one of them (square) is bipartite, 
and the other one (triangular) is not. 

On the square (bipartite) lattice, the correlations are 
power-law^ and the spectrum is gapless: the elemen- 
tary excit ation s are so called resonons with the dispersion 

On the triangular (nonbipartite) lat tice, the correla- 
tions decay expone ntially with distance^ 13 * 18 * 21 !, and the 
spectrum is gappecP^ES with two types of excitations. 
One sector of excitations is of the vison type (topological 
Z2 vortices) , and the other sector is formed by excitations 
local in terms of dimer operators (those excitations can 
be thought of as composed of an even number of vison 
excitations). 

This difference between the physics of the RK states on 
bipartite and non-bipartite lattices can be traced down 
to a larger set of conserved quantities in bipartite dimer 
models. Namely, dimer configurations on bipar tite lat- 
tices may be mapped onto a scalar height fiel d 20 * 22 * 2 ^, 
whose winding numbers provide integer- valued topologi- 
cal invariants (in the dual variables, those integer winding 
numbers may be related to a U(l) gauge symmetry). On 
the other hand, on non-bipartite lattices, only a Z 2 sub- 
group of winding numbers survives, thus breaking the 
U(l) symmetry down to Z 2 . The above difference be- 
tween bipartite and non-bipartite lattices is well known 
and extensively discussed in liter ature, and we refer the 
reader to other studies for details^SES 



III. INTERPOLATION BETWEEN THE RK 
DIMER MODELS ON SQUARE AND 
TRIANGULAR LATTICES 

In our present work, we analyse an interpolation be- 
tween the RK points on the square and triangular quan- 
tum lattices. We introduce this interpolation in such a 
way that the "RK solvability" for the ground state is 
preserved in the interpolating model. Namely, we con- 
struct the interpolation by introducing dimers along one 
of the diagonal directions of the plaquettes of the square 
lattice (Fig. [I]}) . Then by tuning the fugacity of diago- 
nal dimers, one can interpolate between the two limiting 
cases of the square and triangular lattices. 
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FIG. 1: (a) The interpolated model is defined on the square 
lattice with links added along one of the diagonal directions. 
The dashed line connecting plaquettes a and j3 illustrates 
the point-vison construction 0. (b)-(d) The three possi- 
ble dimer-flip processes. The process (b) is contained in the 
"square" part of the Hamiltonian Q; the processes (c) and 
(d) are included in the "diagonal" part jH}. 

The Hamiltonian of the interpolating model is thus 
chosen to be 

H -^square F -f^diag : (3) 

where 

ff 8q uare = E(H I > " I >) (< | " < |) ( 4 ) 

is the RK Hamiltonian ([!]) on the square lattice with 
v = t = 1 (involving kinetic and potential processes for 
square plaquettes depicted in Fig. |TJb) and 

ffdiag = «^ (| Z0 ) - a*I }) (( 0Z I - M 00 l) 

(5) 

is the term involving diagonal dimers. The sum is taken 
over two types of parallelogram plaquettes depicted in 
panels c and d of Fig. [T] The ground state of this Hamil- 
tonian is then 

i GS > = ^E^ (c)/2 i c >> (6) 

v c 

where D(c) is the number of diagonals in a given config- 
uration c. 

Note that this is a two-parameter interpolation: the 
parameter /i is responsible for the ground state of the 
Hamiltonian, while the parameter k determines the dy- 
namics of diagonal dimers. At fj, = K = 0, the Hamilto- 
nian reduces to the RK point on the square lattice, while 
fj, = k = 1 reproduces the RK point on the isotropic 
triangular lattice. 

The above model was introduced in Ref. where 
a simple variational ansatz for the excitations was pro- 
posed. The properties of the ground state ^ were an- 
alyzed in the earlier work Ref. [T3] From the Pfaffian 
technique^, one easily deduces the ground-state corre- 
lation length £ cx One can expect that at any 
nonvanishing /i and k, the system belongs to the Z 2 
phase with excitations continuously connected to those 
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on the isotropic triangular lattice. This scenario has also 
been proposed in Refs. |5] and [24] on the basis of a field- 
theoretical analysis: such an interpolation can be inter- 
preted as a U(l) gauge theory coupled to a charge-two 
scalar field. In the present paper, we complement this 
general analysis with a microscopic construction of vi- 
son excitations in the interpolating model (|3])-([5]) close 
to the critical point (for small deformation parameters /i 
and k). 

In Ref. 26, variational constructions for dimer 
(resonon) and vison excitations have been proposed. For 
the resonon excitations, a quasi-perturbative construc- 
tion gives a gap Adi m oc n 2 n (which is likely to be an ex- 
act asymptotic behavior of the dimer gap in our model) . 
At the same time, an attempt to construct variational 
vison excitations as plane waves of point visons proposed 
in Ref. 26 produced states with a finite variational en- 
ergy that does not tend to zero in the limit of the square 
lattice (/i, k — y 0) . This suggests (and we will confirm it 
in the next section) that point visons are not good vari- 
ational states, but the actual vison eigenstates extends 
over many lattice spacings (in the limit fi, k — > 0). The 
following section is devoted to constructing such states. 



IV. VARIATIONAL CONSTRUCTION OF 
VISON EXCITATIONS 

To construct variational visons, we first define point vi- 
sons and then dress them with local operators to approx- 
imate eigenstates. To define visons, one takes a contour 
T a p connecting two plaquettes a and f3 of the lattice^ 
and considers the intersection-parity operator (Fig. [lP^, 



lower the variational energy, we should spread the vison 
flux from one plaquette to a certain region of the lattice. 
On the square lattice, the point vison may be written in 
terms of the height field as 



V a = exp(iKh a ) 



(10) 



where K is some constant depending on the normaliza- 
tion convention for the height field h a . The spreading 
of the vison flux can then be achieved by the dressing 
operator 



D = exp 







h a ) 



(11) 



with some weights fp determining the flux redistribu- 
tion. We will further use another expression equivalent 



to Eq. (11 1 on the square lattice, but also directly gener- 



alizablc to the interpolating model (where the height field 
h a is no longer defined). Namely, we use the following 
variational ansatz: 



D = exp 



jklT-jk 



(12) 



where the sum is taken over all lattice links and njk is the 
dimer density on the link (jk) (equal to 1 or in the pres- 
ence/absence of a dimer, respectively). Thus constructed 
variational state has a gauge redundancy: a transforma- 
tion 



A jk i ^ A jk + tpj + ipk 



(13) 



y ^ ^ number of dimcrs intersecting T a p 

does not change the state up to an overall phase 



One easily verifies (i) that this operator is independent 
of the contour T a p connecting two given points, up to 
an overall controlled change of sign; (ii) that the com- 
mutator of V a p with any local dimer Hamiltonian is 
concentrated near the points a and j3; and (hi) that 
V a pVp 7 = V ai . Therefore one can decompose the oper- 
ator V a p into the product of two vison operators, which 
have the structure of Z 2 vortices^, 



Therefore, without loss of generality, we may put Ajk = 
on all diagonal links and thereby fix the gauge. Finally, 
it will be convenient to convert Ajk (defined now on the 
square lattice) into a vector potential Ajk by 



A jk = (-I)* A. 



jk 



(14) 



v aP = v a v p . 



(8) 



The variational energy of the "square" part of the 
Hamiltonian Q in the state Q, ( |T2"| ) is given by 

Square = (V| ff square | V) = P,g ^(1 - COS <j>p) (15) 



Alternatively, one can understand the operator V a as the 
parity operator ^ with the point fj sent to infinity. 

Such a point-vison operator produces a state orthogo- 
nal to the original state, but not an eigenstate. To form 
an eigenstate (or, more precisely, a wave packet of lowest- 
energy eigenstates), one needs to dress a point vison, 



|V) =DV a \G$) 



(9) 



where is the probability of a given square plaquette 
to be flippable (at the RK point on the square lattice, 
= 1/4) 26 and cj>p is the flux at the plaquette /?. This 
flux is composed of the flux of the vector potential Ajk 
and of the flux of the original point vison V a in Eq. j9| : 



= A, 



A 2 ?, + A34 + A A1 + ttS, 



aP - 



(16) 



where D is some operator local in terms of dimers. 

We may suggest (and this suggestion is further con- 
firmed by a variational calculation) that, in order to 



where the numerical indices 1 to 4 label the four sites 
of the square plaquette /3 in the cyclic order. Since we 
assume that the operator D is local in terms of dimers 
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and concentrated in a finite region near the plaquette a, 
the total flux of the vector potential Ajk must vanish. 
Therefore, we find the constraint for the total flux <bn. 



^ = 77 ■ 





(17) 



0), the total vari- 
Its minimization 



On the square lattice (at fj, = k = 
ational energy is given by Eq. (151. 
with the constraint ( fi"7| ) yields the spreading of the flux 
over the whole lattice. Indeed, spreading the flux over 
N plaquettes gives N fluxes 6p cx 1/N, which results in 
the total energy E oc 1/N. This variational construction 
shows that, indeed, in the limit of the square lattice the 
vison excitation becomes unstable and that close to the 
critical point the vison size must be large (much larger 
than the lattice spacing). In this limit, using the fact that 
the optimal fluxes 6p are small, we can expand the cosine 
in Eq. ( 15 ) and write the long- wavelength expression for 



the energy, 



E s , 



P (i) 



(V x A) 2 d 2 x. 



with the total screened flux 



A d x = 7T 



(18) 



(19) 



[the latter integral is taken over a small circle surrounding 
the position x a of the point vison in Eq. ([9| ; this point 
is excluded from the integral (flS])]. 

In the interpolating model (at finite \x and k), the sec- 
ond term in the Hamiltonian ^ stabilizes the vison. In- 
deed, its contribution for our variational state Q, (12) 
can be written, in the long- wavelength limit, as 



Er\\c,& 



A 2 d 2 



Here we assume that the vector potential A is small and 
slowly varying on the length scale of one lattice constant. 
This is a valid assumption far from the vison center. The 
central region gives a contribution of the order Kfj, 2 to the 
vison energy, which can be neglected to the leading order, 
in comparison with the logarithmically larger contribu- 

(2) 

tion from large distances. The coefficient P^' is given by 
the probability to find two parallel non-diagonal dimers 
on a "skew" (parallelogram) plaquette. To the leading 
order in /i, it can be approximated by its value at the 
RK point on the square lattice, P^ 0.0454JPI The 
total variational energy of the vison can now be written 
as 



E — E s , 



P 



(1) 



flip 

2~ 



(V x A) 2 + 



A" 



(21) 



where the new large length scale ("size of vison") 
emerges, 



P 



(i) 



\ 8P< 2 



(22) 



iiip 



Remarkably, this variational problem is equivalent to 
that of a vortex in a type-II superconductor (wi th L cor- 
responding to the London penetration lengthJpSESI j n 
terms of the "magnetic field" defined as 



B = V x A, 



(23) 



the variational problem (21 1 is equivalent to minimizing 
the energy 



E = 2 / 
with the total-flux constraint 

Bd 2 X = 7T 



B 2 + L 2 V x B 



(24) 



(25) 



One can find a centrally-symmetric solution to this 
variational problem, 



B(R) 



(26) 



where K is the modified Bessel (Macdonald) function 
and R is the distance from the vison center. The energy 
corresponding to this solution is 

E=-^\nL = -nP^w 2 ln(^ 2 ) . (27) 



Of course, this result for the energy, as well as the long- 
wavelength approximation (18) - (21) and 



(20) 

is only justified if L ^> 1. The expression 



24) - (261, 



26 1 for the 



variational solution is also only valid for R^> 1. 



V. DISCUSSION OF RESULTS 

In the preceding section, we have constructed a varia- 
tional vison excitation in a near-critical RK-type dimcr 
model interpolating between the RK points on the square 
and triangular lattices. While, rigorously speaking, our 
construction provides an upper bound for the gap in the 
vison excitation sector, we believe that it captures the 
main qualitative feature of vison excitations (spreading 
the vison flux over a large region of lattice) , and therefore 
conjecture that it also gives the exact functional depen- 
dence on the deformation parameters /i and n. Below 
we discuss several interesting aspects and implications of 
our result. 

First of all, our variational construction is only based 
on the large length scale L 3> 1, which is equivalent to 
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FIG. 2: (a) The two elementary translations of the lattice, 
which anticommute for vison states, (b) The k space for vison 
excitations. The labeling of the points A, B, and C is the 
same as in the isotropic triangular lattice case considered in 
Ref. I14I (see Fig. 3 there). In our interpolating model, points 
C\ and C2 are not degenerate: they only become degenerate 
in the triangular lattice limit (k — ^ = 1). Close to the critical 
point (square lattice), the vison minima are located at points 
Ci . Note that the vison Brillouin zone is only one half of the 
usual dimer Brillouin zone, and therefore it contains two Ci 
points. 



Kfi <C 1. Thus we do not, in fact, require the individual 
smallness of the parameters k and /i, but only of their 
above combination. 

Second, we find that the (variational) vison gap closes 

as 



-Evison oc — In L 



1 

I? 



(28) 



at the critical point. This behavior should be compared 
to the variational analysis of the dimer gap in the same 
model performed in Ref. 26] There it was found that the 
dimer gap closes as 



1 

1? 



(29) 



i.e., it is always smaller close to the critical point. This 
suggests that visons tend to form bound pairs thus can- 
celing the logarithmic term in the energy (see also the 
discussion of two vison signs below) . 

Third, our variational vison can have two different 
signs, which correspond to the two opposite signs of the 
smeared flux B. While there is only one species of the 
point vison (fluxes 7r and — tt on one plaquette are identi- 
cal), there is a freedom of smearing this flux as either pos- 
itive or negative over many plaquettes (of course, there 
are also options to smear fluxes ±37t, ±57t, etc., but they 
are obviously higher in energy). As a result, we can con- 
struct two variational visons differing by the sign of the 
field B. This result may appear surprising in view of 
our expectation of the Z2 structure of vison excitations. 
This paradox can be resolved by analysing the properties 
of our variational visons with respect to the translational 
symmetries. 



As follows from their construction, visons live on a frus- 
trated lattice^S'. In terms of translational symmetries, 
this implies that the elementary translations in the x\ 
and x 2 directions (along the links of the square lattice, 
Fig. [2^) anticommute: 



TiT 2 



-T 2 T X 



(30) 



Therefore, to define a wave vector for vison eigenstates, 
one needs to double the unit cell. In terms of k vec- 
tors, this implies the period of tt in the directions k\ 
and k 2 for the dispersion relation. More precisely, to 
the four k points (fci, k 2 ), (ki + tt, k 2 ), (hi, k 2 + tt), and 
(fci+7T, k 2 +Tr), there correspond two linearly independent 
vison eigenstates degenerate in energy. 

By analyzing the action of the translation operators 
on our variational vison states, we conclude that these 
states belong to the vicinity of the wave vectors gener- 
ated by k — (0,0). Thus the two vison states of op- 
posite fluxes correspond to the two degenerate minima 
of the vison dispersion relation. Extending the notation 
of Ref. HU these minima can be denoted as C\ (which 
become degenerate with points C 2 in the isotropic tri- 
angular limit, see Fig. ^p). Note that these points are 
different from the energy minima on the isotropic trian- 
gular lattice (k = fj, = 1), which are located at points 

Remarkably, the existence of visons of opposite signs 
can qualitatively explain the formation of vison bound 
states (dimer-like excitations) found numerically in 
Ref. 1141 Indeed, superimposing visons of opposite signs 
at a short distance (much shorter than L) compensates 
the long-distance logarithmic contribution to the energy 
(28), thus producing the energy of the bound state (29). 



Fourth, an analogy with type-II superconductors can 
be drawn not only in the variational form of the energy 
(21) or (24), but also in the existence of two different 
length scales. The vison size L in our problem obvi- 
ously corresponds to the London penetration length in 
the theory of superconductivity. On the other hand, 
there is a second length scale £ (the ground-state cor- 
relation length), which resembles the superconducting 
coherence length. In superconductors, the interplay of 
the two length scales leads to important physical conse- 
quences (difference between type-I and type-II supercon- 
ductors). We may therefore conjecture that in quantum 
dimer models some novel effects depending on the rela- 
tion between the two lengths may also be possible. Note, 
however, that this analogy is not complete: in our vison 
construction, the length scale £ does not enter, and only 
the length L affects the variational state and its energy. 
In fact, the lower cut-off of the logarithm in the vison 



energy ( 28 ) is given by the lattice spacing rather then 
£. It would be interesting to explore further improve- 
ments of our variational ansatz, which may bring in the 
length scale £ (by analogy to the core of superconducting 
vortices). We leave this interesting question for future 
study. 

Another interesting question that remains unexplored 



(> 



is the dynamics of visons in our model. While exact 
vison eigenstates must carry a well-defined wave vector 
k, our variational state is localized in space and thus 
corresponds to a wave packet of width Afc ~ 1/L <C 
1 in the reciprocal space. We believe that its energy 
accurately represents the bottom of the vison band, but 
the question of the vison mass remains unresolved in our 
construction. 

While the study reported in the present paper relates 
to a very specific model and to a non-generic phase tran- 
sition between a RVB phase and a U(l) critical point, 
we believe that some of our findings may be generalized 



to other systems and situations. In particular, the ex- 
istence of two different length scales (£ and L) is likely 
to be a general property of the RVB phase. The diver- 
gence of the vison size may also occur at other types of 
second-order phase transitions. We therefore hope that 
the example considered in the present paper will provide 
a guidance for future studies of the RVB phase in various 
systems and models. 

We are grateful to G. Blatter, M. Feigelman, R. Moess- 
ner, D. Poilblanc, and D. Sadri for helpful discussions and 
comments on this work. 
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